Absence of low temperature anomaly in the Debye- Waller factor of solid "^He down to 
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The mean square atomic displacement in hep-phase solid ''He has been measured in crystals with 
a molar volume of 21.3 cm"^. It is temperature independent from 1 K to 140 mK, with no evidence 
for an anomaly in the vicinity of the proposed supersolid transition. The mean square displacement 
is different for in-plane motions (0.122 ± 0.001 A^) and out-of-plane motions (0.150 ± 0.001 A^). 

PACS numbers: 67.80.-s, 67.80.Cx, 61.12.Ld 



I. INTRODUCTION 
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The recent experiments by Kim and Ghan^ on hexag- 
onal close packed (hep) solid *He have re-opened inter- 
est in the old question of supersolidity. In the 1970s, 
several theorists considered whether the formation of a 
Bose-Einstein condensate was possible in the presence of 
a crystalline lattice^i^i^ and Leggetij^ predicted a signa- 
ture in the rotational inertia. Kim and Chan observed 
just such a shift, and this has subsequently been repro- 
duced by several groups (e.g. Ref. d), although the origin 
and mechanism of this effect remain under discussion. 
Experimental estimates of the onset of this phenomenon 
range from between ~ 80 mK to 800 mK, dependent 
on the percentage of ^He presentj^ An acoustic anomaly 
was observed by Ho et alX in the acoustic attenuation 
in solid ''He crystals, close to 200 mK. This was ascribed 
to a change in the behaviour of defects in the crystal, 
and one school of thought holds that defects may play an 
important role in supersolid formation.™ 

To provide information on the the bulk behavior of 
solid ^Hc in this temperature range, we have carried out 
neutron diffraction experiments on solid "^He crystals to 
look for traces of a supersolid transition in the bulk mate- 
rial. The neutron diffraction cross section is proportional 
to the square of the Fourier transform of the nuclear den- 
sity distribution in space and time. If the nuclei form a 
crystal lattice, then this Fourier transform yields Bragg 
peaks with an intensity that is proportional to the unit 
cell structure factor and the Debye- Waller factor, from 
which the mean square displacement of the ^He nuclei can 
be extracted. Helium forms a quantum crystal and so the 
mean square displacement is greater than the Lindemann 
criterion due to large zero-point fluctuations Depending 
on the model used to describe supersolidity, this quantity 
is liable to change in the supersolid state. As an aside, 
we note that even non-solid samples may have a periodic 
density distribution and hence Bragg peaks, and in this 
case, the effect on the mean square displacement would 



II. EXPERIMENTAL DETAILS 

The He crystals were prepared in a cylindrical stain- 
less steel sample chamber (height 32.1 mm, diameter 33.5 
mm) equipped with quartz transducers to monitor ultra- 
sound propagation through the growing crystal. ^He gas 
was supplied and pressurized through a stainless steel 
capillary located at the top of the sample chamber. A 
mix containing a slightly elevated amountii of ^He (40 
ppm) was used. From Ref. [H this should place the start 
of the transition temperature at ~ 550 mK, although the 
transition reported in Ref. remains close to 200 mK. 

The sample chamber was cooled using an Oxford In- 
struments dilution refrigerator, wherein thermal contact 
to the mixing chamber was established using a Cu wire 
attached to the bottom of the sample chamber. The sam- 
ple chamber supported growth of multiple crystallites 
with different orientation as opposed to a single crys- 
talline sample. 

The dilution refrigerator used for the measurements 
had difficulty cooling to 1 K when the filling capillary to 
the sample chamber was filled with liquid. Therefore, the 
solid had to be formed at higher temperature and pres- 
sure. The sample examined here was formed by starting 
at pressure p = 45 bar and temperature T ~ 2 K. Solid 
''He was then grown using the 'blocked capillary tech- 
nique', in which a solid plug is formed in the capillary as 
the system is cooled so that the molar volume of the sam- 
ple remains constant as it cools along the melting curve. 
This means that our sample passed through the region 
of bee phase on the melting curvCfS. which presumably 
contributed to the formation of strained crystallites. 

The neutron data were collected using the MAPS in- 
strument at the ISIS Facility, Rutherford Appleton Lab- 
oratory, Oxford. MAPS is a time-of-flight (TOE) spec- 
trometer equipped with approximately 16 m^ of posi- 



2 



tion sensitive detectors located 6 m from the sample. 
The instrument was operated in TOF Laue mode with 
a white pulsed incident beam and the sample fixed in 
place. Diffraction patterns were obtained over the tem- 
perature range 140 niK - 800 mK. At the end of the 
experiment the sample can was evacuated and an empty 
can measurement was acquired at 1 K. A standard vana- 
dium sample was also measured to obtain information on 
the wavelength dependent detector efficiency. 



III. DATA ANALYSIS 

To assess the mean square displacement, multiple 
peaks from the same crystallite are required, preferably 
over a large range of momentum transfer. Due to the 
presence of multiple crystallites in the sample can, it was 
impractical to identify peaks from the same crystallite 
at different locations on the detector. However, because 
we used a pulsed white incident beam, higher-order re- 
flections from the same crystallite appeared in the same 
detector pixels, separated by the recorded time-of-flight. 
For these higher-order reflections to be visible, the princi- 
pal reflection had to be sufficiently strong. 45 reflections 
were indexed on the basis of their d-spacings, assuming 
an hep structure. Peaks of the types (002), (100), (101) 
and (102) were readily identified, but only peaks of the 
type (002) and (100) were strong enough for the corre- 
sponding higher-order reflections (004) and (200) to be 
visible at half the TOF. 

The molar volume was determined using the range of 
plane spacings obtained at this point, even though the 
peaks were not from the same crystallites. We found 
a ~ 3.68(1) A and c = 6.03(1) A, corresponding to a 
molar volume of 21.3(1) cm'^. The c/a ratio is 1.638(5) 
as compared with the ideal hep ratio of 1.633. This is 
consistent with expectations from the temperature and 
pressure during the crystal growth process.— 

Once these peaks were identified, the TOF spectra in 
the corresponding pixels were normalized to the incident 
flux on the sample and the empty can measurement sub- 
tracted. At small momentum transfers, where the ma- 
jority of scattering from the can and cryostat was found, 
a small amount of background scattering remained after 
this subtraction. An additional subtraction was made 
when integrating over reflections by taking a slice from 
the same pixel over a d-spacing range close to that of 
interest. At higher momentum transfers, the initial sub- 
traction is sufficient. 

The integrated intensities were then corrected for (i) 
angular effects, (ii) detector efficiency and (iii) the wave- 
length dependence of the incident beam fiux. To estab- 
lish the first correction factor, consider a diffraction peak 
with a scattering angle 29b acquired on a 2D plate a dis- 
tance Z (= 6.0031 m) from the sample. The centre of the 
detector plate coincides with the centre of an orthonor- 
mal coordinate system where x and y are perpendicular 
directions in the plane of the detector, and z is perpendic- 



ular to the plate. The reflection is observed at point (x,y) 
on the plate, and distance from the centre of the sample 
to this point is i? = (x^ + + Z^)!/^. The direction of 
the scattered beam can be expressed as a function of two 
angles, corresponding to the latitude and longitude (la- 
belled ip and (f>). Using these angles, the position of the 
Bragg peak is _R(cos ijj sin 0, — sin -0, cos ip cos 0), keeping 
the (x, y, z) coordinate system. The scattering vector Q 
is (^,^,fc(|r - 1)) and so = 2fc2(l - cos cos 0) = 
4fc^ sin^ 6b, where k is the magnitude of the wavevector. 
The number of neutrons counted in the Bragg peak is 
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where $ is the incident flux on the sample (neutrons / 
unit area / unit time), and Af2 is the solid angle sub- 
tended by the sample. We convert this to an integral 
form; the solid angle can be written as {dS ■'R)/R^ = 
dxdy Z/R^. The incident fiux can be written as a func- 
tion of the incident wavelength - any variation from unity 
here is accounted for when the vanadium scans are con- 
sidered. 



$ = /(A)dA = 



A^/(A) 
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dk 



The neutron scattering cross-section is 



r) 
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where N is the number of unit cells divided by the vol- 
ume of a unit cell, and is particular to each crystallite, 
and |F(t)P is the structure factor of the Bragg peak. 
We therefore arrive at the following expression for the 
intensity of the Bragg peak at reciprocal lattice vector t. 

= (27r)37V|F(r)p^!M^y' j j 5{Q-t) dxdydk 

, (4) 

Integrating over momentum space, the Bragg peak inten- 
sity is given by 



lB=2N\F{T)\''dl^i\^I{X). 



(5) 



This expression applies equally to all of the detector 
banks. The value of dhki for the main peaks was obtained 
by fitting the peak profile as a Gaussian and taking the 
central value. This value was then halved for the relevant 
higher-order peak. To correct for the other factors, the 
vanadium scan was used. 

Vanadium is a uniform scatterer, and so 

/vana = J^/(A) AAAArf -J 

= ^^^cos^26'Bsin6'B/(A)Ad (6) 

where = 5.08-10"^* m^ is the incoherent cross-section 
for vanadium^ii Z = 6.0031 m, and AAd = iVpAxA?/ is 
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Figure 1: The mean square displacement, (u^), of ''He licp 
crystal, as measured for the three selected peaks as a function 
of temperature. 



Figure 2: The crystallite dependent quantity In A^, as mea- 
sured from *He hep crystal for the three selected peaks as a 
function of temperature. 



the relevant area of the detector {Np is the number of 
pixels). For one pixel Ax = 0.025 m (detector tube diam- 
eter) and Ay = 0.015 m (vertically resolvable distance). 
Eq. [S] allows /(A) to be extracted from the white beam 
vanadaium diffraction data. 

In the harmonic approximation, the structure factor 
|F(t)P is related to the mean-square displacement as 
follows 

F(t) = bHe cxp(zx • d,) exp{~Wd, (t)). (7) 

i 

Here is the position of the i*^ atom in the unit cell, 
bne = 3.26-10~^'^ m is the scattering length of '^He^ii 
Wdi{T) ~ {{t ■ u)^), and u is the displacement from the 
average periodic lattice for a specific nucleus For the 
He hep crystal, the lowest order expansion of the mean 
square displacement can be broken down into in-plane 
and out-of-plane (|| to c) components, as follows 

((x . nf) = AiuDir ■ cf) + {ul){l [t ■ cf))]. (8) 

Previous studies have found this expression to be quite 
adequate for describing the observed dataj ^'^i^'^ 

Hep He has two atoms per unit cell, and due to sym- 
metry equivalence Wd^ has no site dependence, and so 
can be pulled out of the summation over i. After making 
these corrections, one finds that 

|exp-(W^dJP==^ (9) 

where /corr = C [I b I Ivans.) IS the corrected intensity of 
each peak, and C is the total correction factor, obtained 
from Eqs.liElandlll 
We therefore have 

ln/eo„ = In TV - 202[(y|)(Q . cf) + («i)(l - (Q • cf))\ 

(10) 



and so the mean square displacement can now be deter- 
mined from a simple linear fit of intensity vs . 

The level of thermal diffuse scattering, calculated fol- 
lowing the method of Popa and Willisi^ for TOF neutron 
diffractometry, using the phonon velocities measured by 
Minkiewicz et al.,— was found to be negligible. 

IV. RESULTS 

Three peaks met the conditions established above, and 
values for (u^) were obtained at each temperature consid- 
ered. Figures [T] and show the temperature dependence 
of (u^) and In TV. The data are tabulated in the appendix. 
No thermal anomaly indicative of a phase transition is 
observed, consistent with recent path-integral computa- 
tions that indicate that bulk solid He crystals cannot be 
supcrsolidfii Table |T] averages over all of the temperatures 
assuming no temperature dependence, (u^) is higher for 
the (OOL) peaks than for the (HOO) peaks indicating some 
anisotropy. The similar values of In N indicate similar 
crystallite sizes. Consideration of the actual peak posi- 
tions on the detector indicate that the second and third 
peaks are (100) type peaks that could be from the same 
crystallite, while the (002) type peak is from a different 
crystallite. 

Figure [3] shows the temperature dependence of the lat- 
tice parameter a [c] as measured from the (100) [(002)] 
reflection. Although the lattice parameters for a given 
crystallite remained fairly stable as a function of temper- 
ature, due to different annealing conditions and difficul- 
ties with temperature stability between different runs, 
there was some change in the values obtained. This is 
most clearly obvious in the upper panel of Figure[3l where 
two regimes are visible. The upper set of points, which 
differ from the rest by 1 part in 6000, correspond to condi- 
tions where the stress on the crystallite was believed to be 
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Figure 3: The temperature dependence of the c and a lattice 
parameters from two of the crystalhtes studied in the text. 



slightly different to that present for the other measured 
points. These changes had little impact on the observed 
mean square displacement, indicating that we are ob- 
serving the (temperature- and volume-independent) zero- 
point motion here. This is in agreement with the obser- 
vations made at a lower molar volume by Adams et al^ 
from 0.07 - 0.4 K. 



V. DISCUSSION 





Peak 


d-spacing (A) 


203 (°) 




IniV 


1 


002 


3.016 


20.6 


0.150 ± 0.001 


58.560 ± 0.002 


2 


100 


3.188 


18.8 


0.118 ± 0.001 


58.716 ± 0.003 


3 


100 


3.186 


19.7 


0.125 ± 0.001 


58.636 ± 0.002 



Table I: Characteristic properties of the peaks examined on 
MAPS. 



Molar 


Temp. 




Peak 


Reference 


volume (cm'^) 


(K) 


(A^) 


type 




11.01 


15 


0.0593(1) 


mixed 


Ref. J_3 (x-ray) 


12.06 


5.8 


0.0466(3) 


mixed 


Ref. JA (neutron) 


12.12 


14.8 


0.0563(14) 


mixed 


Ref. 18 (x-ray) 


12.13 


14.8 


0.0513(10) 


mixed 


Ref. _18 (x-ray) 


15.72 


5.8 


0.0861(9) 


mixed 


Ref. JA (neutron) 


20.9 


0.7 


0.1537(7) 


(OOL) 


Ref. 29 (x-ray) 


21.3 


<1 


0.150(1) 


(OOL) 


this work (neutron) 


21.3 


<1 


0.122(1) 


(HOO) 


this work (neutron) 



Table II; Published values for {u^) in solid *He obtained by 
x-ray and neutron diffraction. In all cases, the harmonic 
approximation is assumed and for the mixed peak data, no 
distinction between in-plane and out-of-plane displacements 
were made. 



the temperature dropped in the fee solid phase. We see 
no temperature dependence, which implies that we are 
looking at quantum fluctuations. 



Table mi lists pubhshed data for (m^) at different molar 
volumes. The form for the Debye- Waller factor differs 
in some of these references, and so all values have been 
converted to the definition of (u^) given above. These 
values compare well with theoretically calculated num- 
bers; the value for a crystal of molar volume 15.72 cm'' 
is reproduced very well by path-integral Monte Carlo 
techniqueS)^^ and the value at 20.9 cm^ is close to that 
found by Whitlock et al. using the Green's function 
Monte Carlo methodi^ Some estimates were obtained 
from inelastic neutron scattering studies of the phonons, 
but the validity of these estimates have been called into 
doubt by Burns and Isaacs,— who claim that the phonon 
data suffer from significant multiple scattering and arc 
not included. 

We have observed a certain amount of anisotropy be- 
tween in-plane and out-of-plane displacements. The in- 
plane value is 20 % smaller than the out-of-plane value, 
indicating that atomic motion is easier out-of-plane. This 
is not unexpected for an hep structure, particularly as the 
close packing is not perfect, and the qualitative behavior 
agrees with that calculated by Reese et al^ from phonon 
frequencies. The previous x-ray study by Venkatara- 
man and Simmonsi^ did not report any anisotropy in 
the hep phase, but the reflections studied are not noted. 
Venkataraman and Simmons noted a decrease in (u^) as 



The lattice parameter measurements allow us to com- 
ment on the model proposed by Anderson et al^, that 
posits a dependence for the vacancy concentration, as 
opposed to a classical thermally activated model. Fraass 
et al^ used the change in lattice parameter with temper- 
ature as a measure of the vacancy concentration. If the 
actual mass density of He atoms is constant as a function 
of temperature, this approach is valid, and the caveats as- 
sociated with this analysis are outlined in Ref. [l^- From 
Fraass et al. 's data, the two models could not be sep- 
arated, but with the extension to temperatures below 
0.8 K provided here, the thermally activated model is 
favoured (see Fig. 3]) , and fitting the complete data set 
gives an formation energy of (8.6 ± 0.1) K. We note that 
the samples compared here had slightly different molar 
volumes, and that in our sample chamber, there were sev- 
eral crystallites and so the volume of a given crystallite 
may have varied as a function of temperature. 

As a complement to this, we note that acoustic at- 
tenuation data collected by one of us^^ provide an ad- 
ditional confirmation of this exponential dependence, as 
the acoustic attenuation is determined by the combina- 
tion of phonon and dclocalizcd vacancy contributions. 
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VI. CONCLUSIONS 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Temperature (K) 



Figure 4: The temperature dependence of the a lattice pa- 
rameter as measured here (closed points) and as calculated 
from Fig. 2 of Ref. [S^ (open points) , assuming that the tem- 
perature independent part of a is the same for both data sets. 
The lines are based on the different models for vacancy con- 
centration: the solid line is the exponential from the thermal 
activation model, with an activation energy of 8.6 K; the bro- 
ken lines are fits to the model proposed in Ref. [22 - a 
dependence. The dotted line is a fit to the data from Ref. [S^ 
only and the dashed line includes the data measured in this 
paper. 



The principal result of this paper is that the mean 
square displacement does not change over the tempera- 
ture range measured here. Within the precision of these 
data, there is no apparent temperature dependence at all 
over the temperature range studied, i.e. there is no indi- 
cation that the supersolid transition, if it exists, affects 
the crystalline lattice or zero-point fluctuations. We have 
also observed a measure of anisotropy between in-plane 
and out-of-plane motions. 
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Appendix: DATA TABLES 
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OOL HOOa HOOb 



Temp. 


{ 




error 


{ 


u') 


error 




error 


(K) 


( 


k^) 


(A^) 


(A^) 


(A^) 


(A=) 


(A=) 





140 





151 


0.004 





119 


0.007 





134 


0.006 





165 





165 


0.009 





136 


0.015 





130 


0.008 





165 





152 


0.004 





120 


0.006 





130 


0.005 





165 





150 


0.004 





123 


0.007 





126 


0.004 





175 





159 


0.006 





115 


0.008 





131 


0.006 





178 





146 


0.004 





121 


0.009 





127 


0.006 





180 





155 


0.005 





136 


0.012 





129 


0.006 





180 





150 


0.003 





118 


0.006 





125 


0.004 





195 





151 


0.005 





123 


0.009 





116 


0.005 





195 





160 


0.008 





129 


0.014 





121 


0.007 





200 





153 


0.005 





114 


0.007 





133 


0.007 





200 





164 


0.007 





118 


0.008 





127 


0.006 





200 





146 


0.004 





144 


0.015 





143 


0.008 





200 





146 


0.003 





116 


0.006 





126 


0.004 





210 





149 


0.005 





127 


0.010 





130 


0.006 





215 





150 


0.005 





112 


0.007 





132 


0.007 





220 





147 


0.004 





129 


0.010 





122 


0.005 





220 





148 


0.004 





126 


0.010 





127 


0.006 





220 





154 


0.005 





112 


0.007 





137 


0.007 





220 





141 


0.003 





115 


0.006 





121 


0.004 





240 





152 


0.004 





106 


0.005 





123 


0.004 





300 





154 


0.004 





115 


0.006 





121 


0.004 





350 





153 


0.004 





123 


0.007 





125 


0.004 





400 





147 


0.003 





120 


0.006 





122 


0.004 





450 





147 


0.003 





119 


0.006 





115 


0.003 





500 





161 


0.005 





123 


0.007 





132 


0.005 





500 





181 


0.011 





108 


0.006 





122 


0.005 





600 





154 


0.004 





135 


0.009 





116 


0.003 





800 





152 


0.004 





121 


0.006 





132 


0.005 



Table III; The mean square displacement (u^) and associated 
error for all temperatures measured, as measured for three 
selected Bragg peaks from hep crystals of "^He. 
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OOL HOOa HOOb 



Temp. 


InN 


error 


InN 


error 


InN 


error 


(K) 


(A^) 








(A^) 


(A^) 


(A^) 


0.140 


58.574 


0.009 


58 


727 


0.015 


58.721 


0.013 


0.165 


58.834 


0.021 


58 


636 


0.013 


58.623 


0.011 


0.165 


58.560 


0.008 


58 


748 


0.015 


58.649 


0.009 


0.165 


58.577 


0.008 


58 


720 


0.013 


58.679 


0.010 


0.175 


58.689 


0.018 


58 


845 


0.031 


58.673 


0.016 


0.178 


58.552 


0.007 


58 


705 


0.013 


58.636 


0.009 


0.180 


58.524 


0.007 


58 


691 


0.012 


58.645 


0.009 


0.180 


58.476 


0.006 


58 


701 


0.012 


58.601 


0.008 


0.195 


58.564 


0.008 


58 


624 


0.010 


58.616 


0.009 


0.195 


58.583 


0.008 


58 


701 


0.012 


58.606 


0.008 


0.200 


58.578 


0.008 


58 


758 


0.014 


58.638 


0.009 


0.200 


58.520 


0.007 


58 


733 


0.013 


58.607 


0.009 


0.200 


58.518 


0.007 


58 


728 


0.013 


58.559 


0.008 


0.200 


58.644 


0.010 


58 


760 


0.014 


58.692 


0.011 


0.210 


58.587 


0.008 


58 


857 


0.018 


58.564 


0.008 


0.215 


58.570 


0.008 


58 


742 


0.014 


58.686 


0.010 


0.220 


58.532 


0.009 


58 


910 


0.030 


58.782 


0.017 


0.220 


58.688 


0.014 


58 


712 


0.017 


58.657 


0.012 


0.220 


58.555 


0.009 


58 


781 


0.020 


58.680 


0.013 


0.220 


58.560 


0.010 


58 


664 


0.015 


58.700 


0.014 


0.240 


58.543 


0.009 


58 


772 


0.020 


58.661 


0.012 


0.300 


58.538 


0.009 


58 


791 


0.021 


58.615 


0.011 


0.350 


58.598 


0.011 


58 


657 


0.015 


58.731 


0.015 


0.400 


58.585 


0.010 


58 


672 


0.015 


58.703 


0.014 


0.450 


58.642 


0.012 


58 


684 


0.016 


58.690 


0.013 


0.500 


58.527 


0.009 


58 


728 


0.018 


58.655 


0.012 


0.500 


58.604 


0.011 


58 


849 


0.025 


58.673 


0.013 


0.600 


58.569 


0.010 


58 


739 


0.019 


58.574 


0.010 


0.800 


58.646 


0.016 


58 


788 


0.028 


58.606 


0.014 



Table IV: The crystallite dependent quantity In N and asso- 
ciated error for all temperatures measured, as obtained from 
''He hep crystals for three selected peaks. 



